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Very Important 

This article is accompanied by the Maple packages 

• http: //www. math. rutgers . edu/~zeilberg/tokhniot/RITSUF , 

• http: //www. math. rutgers . edu/~zeilberg/tokhniot/RITSUFwt , 

• http: //www. math. rutgers . edu/~zeilberg/tokhniot/ARGF , 

to be described below. In fact, more accurately, this article accompanies these packages, written by 
DZ and the many output files, discovering and proving deep enumeration theorems, done by SBE, 
that are linked to from the webpage of this article 

http : / / www . math . rutgers . edu/~zeilberg/mamarim/mamar imhtml/r it suf . html 
How It All Started: April 5, 2012 

During one of the Rutgers University Experimental Mathematics Seminar dinners, the name of Don 
Knuth came up, and two of the participants, David Nacin, who was on sabbatical from William 
Patterson University, and first-year graduate student Patrick Devlin, mentioned that they recently 
solved a problem that Knuth proposed in Mathematics Magazine{Kn]. The problem was: 

1868. Proposed by Donald E. Knuth, Stanford University, Stanford, CA. 

Let n > 2 be an integer. Remove the central (n — 2) 2 squares from an (n + 2) x (n + 2) array of 
squares. In how many ways can the remaining squares be covered with An dominoes? 

As remarked in the published solution in Math. Magazine, the problem was already solved in the 
literature by Roberto Tauraso[T]. The answer turned out to be very elegant: 4(2i^ + (— l)") 2 . 

David and Patrick, as well as the solution published in Math. Magazine, used human ingenuity 
and mathematical deduction. But, as already preached in [Zl] and [Z2], the following is a fully 
rigorous proof: 

"By direct counting of tilings, the first 16 terms (starting at n = 2) of the enumerating sequence 
are 
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47155859716, 323212716324, 2215328606404, 15184099435684, 104073336269956 . 

But so are the first 16 terms of the sequence {4(2F% + (— l) n ) 2 }. Since the statement is true for 
the first 16 terms, it must be true for all n>2. QED!" 

In order to justify this "empirical" proof, all we need to say is that both sides are obviously C-finite 
sequences whose minimal recurrences have order < 8, hence their difference is a C-finite sequence of 
order < 16 and hence if it vanishes for the first 16 terms, it always vanishes. But first let's remind 
ourselves what are C-finite sequences. 

C-finite Sequences 

Recall that a C-finite sequence {a(n)}™ =0 is a sequence that satisfies a homogeneous linear- 
recurrence equation with constant coefficients. It is known (but not as well-known as it 
should be!) and easy to see (e.g. [Z2],[KP]) that the set of C-finite sequences is an algebra. 
Even though a C-finite sequence is an "infinite" sequence, it is in fact, like everything else 
in mathematics (and elsewhere!) a finite object. An order-L C-finite sequence {a(n)}^L ^ s 
completely specified by the coefficients ci, C2, . . . , cl of the recurrence 

a(n) = c\a{n — 1) + C2<a(n — 2) + . . . + C£,a(n — L) , 

and the initial conditions 

o(0) = d\ , ... , a(L — 1) = d,L 
So a C-finite sequence can be coded in terms of the 2L "bits" of information 

[[d 1 ,...,d L ],[c 1 ,...,c L \] . 
For example, the Fibonacci sequence is written: 

[[0,1], [1,1]]. 

Since this ansatz (see [Zl]) is fully decidable, it is possible to decide equality, and evaluate ab 
initio, wide classes of sums, and things are easier than the holonomic ansatz (see, e.g., [Z2]). The 
wonderful new book by Manuel Kauers and Peter PaulefKP] also presents a convincing case. 

Rational Generating Functions 

Equivalently, a C-finite sequence is a sequence {a(n)} whose ordinary generating function, 

Yl^Lo a ( n ) zn , is rational and where the degree of the denominator is more than the degree of the 
numerator. These come up a lot in combinatorics and elsewhere (e.g. formal languages). See the 
old testament[St], chapter 4, and the new testament[KP], chapter 4. 
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Why is the Number of tilings of the Knuth Square-Ring C-finite? 

Each of the four "corners" (2x2 squares) is tiled in a certain way, where either the tiles covering it 
only cover it, or some tiles also cover neighboring cells not in the corner-square. There are finitely 
many such scenarios for each corner-square, hence for the Cartesian product of these scenarios. 
Having decided how to cover these four corner squares, one must decide how to tile the remaining 
four sides, each of which is either a 2 by n rectangle, or one with a few bites taken from one or both 
ends. By the transfer matrix method ([St], 4.7), it follows a priori that each of these enumerating 
sequences is C-finite, hence their product, and hence adding up these finitely many sequences. It is 
also easy to see that a rigorous upper bound for the order is 8, so the above "empirical" approach 
is justified. 

But we still need to be able to compute the first 16 terms. If we had a very large and very fast 
computer, we can actually construct all the tilings, and then count them, but 104073336269956 is 
a pretty big number, so we need a more efficient way. 

A More Efficient Way 

Suppose that you are given a set of unit-squares (let's call them cells) and you want your computer 
to find the number of ways of tiling it with dominoes. You pick any cell (for the sake of convenience 
the left-most bottom-most cell), and look at all the ways in which to cover it with a domino piece. 
For each of these ways, removing that tile leaves a smaller region, thereby getting an obvious 
dynamical programming recurrence. 

Procedure NT(R) in the Maple package RITSUF computes the number of domino tilings of any set of cells. 

Using this method, the first-named author solved Knuth's problem in 15 seconds, by typing (in a 
Maple session, in a directory where RITSUF has been downloaded to) 

read RITSUF: SeqFrameCsqDirect(2,2,2,2,20,t) ; . 

We will soon try to solve analogous problems for fatter frames, but then we would need to be even 
more efficient, and using this more efficient method, the same calculation would take less than half 
a second, typing: 

read RITSUF: SeqFrameCsq(2,2,2,2,20,t) ; . 
An even More Efficient Way 

Our general problem is to find an efficient automatic way to compute the C-finite description, 
and/or its generating function (that must always be a rational function) for the number of domino 
tiling of the region, that we denote, in RITSUF, by 

Frame (al ,a2 ,bl ,b2 ,n,n) , 
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that consists of an (al + n + a2) x (61 + n + 62) rectangle with the middle n x n square removed. 



Before describing the algorithm, let us mention that this is accomplished by procedure 

SeqFrameC(al,a2,bl,b2,N,t) 

For example, SeqFrameC (1,3, 3,1, 30, t); yields: 

-9 + 8t + 29t 2 - 10t 3 -7t 4 + 2t 5 



-4- 



(t 2 - 4i + 1) (t 4 - 4t 2 + 1) 



The fifth argument, N is a parameter for "guessing" the C-finite description, indicating how many 
data points to gather before one tries to guess the C-finite description. It is easy to find a priori 
upper bounds, but it is more fun to let the user take a guess, and increasing it, if necessary. 

Mihai Ciucu's Amazing Theorem 

The sequence of positive integers, demanded by Knuth, enumerating the domino tilings of Frame(2, 2, 2, 2, n, n) 
turned out to be all perfect squares. This is not a coincidence! A beautiful theorem of Mihai 
Ciucu[C] tells us that whenever there is a reflective symmetry, the sequence enumerating the dimer 
counting is either a perfect square or twice a perfect square. Since we know that (and even if we 
didn't, we could have discovered it empirically for each case that we are computing), we have to 
gather far less data points. For example, according to the first-named author of this article, the 
number of domino tilings of Frarae(3, 3, 3, 3, n, n) is 2B(n) 2 , where 

^ -29 + 19 1 + 102 t 2 - 32 1 3 - 25 t 4 + 7t 5 

^ B{n)t ~ (t 2 -4t + l)(t 4 -4t 2 + l) 

while the number of domino tilings of -Frame(4, 4, 4, 4, n, ri) is C(n) 2 where 

P(t) 



n=0 

where 



Q(t) 



P(t) = -4(901+2517 1-17574 1 2 -46322 1 3 +112903 1 4 +291045 1 5 -269376 1 6 -741508 1 7 +215233 1 8 +786069 1 9 - 

21836 1 10 - 352896 1 11 - 24137 1 12 + 67487 1 13 + 5874 1 14 - 5056 1 15 - 359 1 16 + 97 1 11 ) , and 
Q{t) = (t- 1) (t + 1) (t 4 + t 3 - 5t 2 + t + 1) (i 4 - lit 3 + 25 1 2 - lit + 1) • 
(t 4 + 7t 3 + 13t 2 + 7t + l) (t 4 -t 3 -5t 2 -t + l) . 



If you want to see the analogous expressions for Frame(5, 5, 5, 5, n, n) and Frame(6, 6, 6, 6, n, n), 
then you are welcome to look at the output file 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/oRITSUF2, 
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that you can generate yourself, ab initio by running (once you uploaded RITSUF onto a Maple 
session) 

http : //www . math . rutgers . edu/~zeilberg/tokhniot/ inRITSUF2 

Our method is pure guessing, but in order to guess, we need to efficiently generate sufficiently many 
terms of the counting sequence. We must start with rectangles of fixed width. 

The Number of Domino Tilings of a Rectangle of a Fixed Width 

Let m be a fixed positive integer. We are interested in a C-finite description, as a function of n, of 
the sequence A m (n), the number of domino tilings of an m x n rectangle. In fact, for this specific 
problem there is an "explicit" solution, famously found by Kasteleyn[Ka] and Fisher & Temperly 
[FT], but their solution only applies to domino tiling, and we want to illustrate the general method. 

Also the general approach, using the transfer-matrix method, is not new as such (see, e.g. [St]), 
but since we need it for counting more elaborate things, let's review it. 

Consider the task of tiling the n columns of an m x n rectangle. Let's label the cells of a given 
column from bottom to top by {1, . . . , m}. When we start, all the m cells of the leftmost column 
are available, so we start with the state {1, . . . , m}. As we keep on going, not all cells of the current 
column are available, since some of them have already been tiled by the previous column. The other 
extreme is the empty set, where nothing is available, and one must go immediately to the next 
column, where now everything is available, i.e. the only follower of is the universal set {1, . . . , m}. 

Let's take m = 4 and see what states may follow the state {1,2,3,4}. We may have two vertical 
tiles {1, 2} and {3, 4}, leaving all the cells of the next column available, yielding the state {1, 2, 3, 4}. 
We may decide instead to only use horizontal tiles, leaving nothing available for the next column, 
resulting in the state 0. If we decide to have one vertical tile in the current column, then 

If it is {1, 2} then both 3 and 4 are parts of horizontal tiles that go to the next column, leaving the 
set of cells {1,2} available. 

If it is {2, 3} then both 1 and 4 are parts of horizontal tiles that go to the next column, leaving the 
set of cells {2, 3} available. 

If it is {3, 4} then both 1 and 2 are parts of horizontal tiles that go to the next column, leaving the 
set of cells {3, 4} available. 

It follows that the "followers" of the state {1,2,3,4} are the five states 

0, {1,2,3,4}, {1,2}, {2,3}, {3,4} . 

Who can follow the state {1,4}? Since only cell 1 and cell 4 are available there can't be any vertical 
tiles, and both must be parts of horizontal tiles, occupying {1,4} of the next column, and leaving 
{2,3} available, so the state {1,4} has only one follower, the state {2,3}. 



5 



Check out procedure Followers (S ,m) in Ritsfuf 



This way we can view any tiling of an m x n rectangle as a walk of length n in a directed graph whose 
vertices are labeled by subsets of {1, ... , m}. The transition matrix of this graph is the 2 m x 2 m 
matrix whose rows and columns naturally correspond to subsets (in RITSUF we made the natural 
convention that the indices of the rows and columns correspond to the binary representations of the 
subsets plus 1). Let's call SN(i) the set of positive integers corresponding to the positive integer 
i. For example, SN(1) is the empty set, SN (10) is {1,4} etc. 

Check out procedure TM(m) in Ritsfuf for the transition matrix. 

Calling this transfer matrix A m , the number of domino tilings of an m x n rectangle is the (2 m , 2 m ) 
entry of the matrix A 7 ^, since we have to completely tile it, the starting state must be {1, ... , m}, 
of course, but so is the ending state, since every thing must be covered, leaving the next column 
completely available. 

Check out procedure SeqRect(m,N) in Ritsfuf for the counting sequence for the number of domino tilings of an Tfl 
by n rectangle for n = 0, 1, . . . , N . 

But not just the (2 m ,2 m ) entry of A^ is informative. Each and every one of the (2 m ) 2 entries 
contains information! A typical (i,j) entry of A^ tells you the number of ways of tiling an m x n 
rectangle where the leftmost column only has the cells in SN(i) available for use, while the rightmost 
column has some tiles that stick out, leaving available for the (n + l)-th column the cells of SN(j). 

Counting the Number of Domino Tilings of a Holey Rectangle 

We want to figure out how to use matrix-multiplication to determine the number of tilings of the 
region 

Frame(ai, a 2 , 61, b2,m, n) 
that consists of an (a\ + m + a 2 ) x (61 + n + b 2 ) rectangle with the middle mxn rectangle removed. 

There are four corner rectangles in our frame: 

• the left-bottom (SW) corner consisting of an a\ x b\ rectangle , 

• the right-bottom (SE) corner consisting of an a\ x b 2 rectangle , 

• the left-top (NW) corner consisting of an a 2 x bi rectangle , 

• the right-top (NE) corner consisting of an a 2 x b 2 rectangle 

If we look at a typical tiling of the region Frame(ai,a 2 ,bi,b 2 ,m,n), and focus on the induced 
tilings of the four corner-rectangles, we get a tiling with (usually) some of the tiles intersecting the 
adjacent non-corner rectangles. 

Indeed, for the left-bottom a\ x 61 corner-rectangle (usually) some of the tiles covering its very top 
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row intersect the very bottom row of the left (West) mxbi rectangle, and (usually) some of the tiles 
covering its very right column intersect the leftmost column of the bottom (South) a\ x n rectangle, 
and similarly for the other three corner rectangles. So suppose that the East m x b\ rectangle has 
already been tiled, with some of its bottom tiles overlapping with our above-mentioned left-bottom 
a\ x b\ corner-rectangle, leaving only some of the cells in the top row available, and after we 
complete the tiling of that left-bottom ai x bi corner-rectangle, we may use-up some of the cells of 
the left column of the bottom (South) ai x n rectangle, and the complement of that occupied set 
is only available for tiling. 

This leads naturally to a transfer matrix between the "states" of one side of a corner-rectangle to 
the states of another side of that corner-rectangle. So let's define RTM{a,b,Si,S 2 ), for positive 
integers a and b, and distinct S\, S 2 € {1, 2, 3, 4} where we make the convention 

l=Up Side , 2=Left Side , 3=Down Side , 4=Right Side , 

that tells you the number of ways of tiling the a x b rectangle where the two sides that are not in 
{Si, S 2 } are "smooth" (i.e. nothing sticks out) and the two sides Si, S 2 may (and usually do) have 
some of their tiles "sticking out" . 

Going counterclockwise starting at the left-bottom (SW) corner, we need to find 

• For the left-bottom (SW) corner consisting of an ai x b\ rectangle we need RTM{ai, b\, 1, 4) , 

• For the right-bottom (SE) corner consisting of an ai x b 2 rectangle we need RTM(a\, b 2 , 2, 1) , 

• For the right-top (NE) corner consisting of an a 2 x b 2 rectangle we need RTM(a 2 , b 2 ,3, 2) , 

• For the left-top (NW) corner consisting of an a 2 x bi rectangle we need RTM(a 2 , 61, 4, 3) 

Like the matrices TM(m) that for each needed m, we only compute once and then record it, (using 
option remember), we also only compute RTM(a,b, Si, S 2 ) once for each needed a,b,Si,S 2 and 
remember it. 

But how to compute RTM(a,b, Si, S 2 )? We first construct, literally, all the domino tilings that 
completely cover the cells of the a x b rectangle where nothing sticks out of the sides that are not 
labelled Si or S 2 , but that may (and usually do) stick out from the sides labelled Si and S 2 . Then 
we look at all the pairs of states, and form a matrix whose entry is the number of stuck-out 
tilings of the a x b rectangle where the "stick-out" state of side Si corresponds to the set labelled 
i, and the "stick-out" state of side S 2 corresponds to the set labelled j. 

See procedure RTM(a,b, SI ,S2) for its implementation in RITSUF. 

It follows that, in terms of the matrices TM{m) and RTM(a,b, Si, S 2 ), the quantity of interest, 
the number of tilings of the rectangular picture-frame Frame(ai,a 2 ,bi,b 2 ,m,n), is 

Trace{RTM(ali,bi, 1, A)TM{al) n RTM(ai,b 2 , 2, l)TM(b 2 ) m RTm(a 2 , b 2 ,3, 2)TM{a 2 ) n RTM(a 2 , bi, 4, 3)TM(&i) ; 
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Since matrix power-raising is very fast, and so is matrix-multiplication, we can quickly crank-out 
sufficiently many terms in the enumerating sequence, and since we know a priori that it is C-finite, 
we can guess a C-finite description (and/or a rational generating function), that is proved rigorously 
a posteriori by checking that the order bounds are right. 

The Bivariate Generating Function 

If you are interested in the discrete function of the two discrete variables m and n, for the number 
of domino tilings of Frame(a\, a 2 ,b 1 , b 2 ,m, n), then it is doubly C-finite, meaning that its bivariate 
generating function has the form P(x, y)/(Qi{x)Q 2 (y)), for some polynomials P(x, y), Qi(x), Q 2 (y). 
Using an analogous method for guessing, after cranking-out enough data, we can get these gener- 
ating functions easily, using procedure 

GFf rameDouble (al,a2,bl,b2,x,y,N) 

in RITSUF. For example, if D(m, n) is the number of domino tilings of Frame(2, 2, 2, 2, m, n), then 

oo oo p/ \ 

m=0n=0 W\X,y) 

P{x,y) = 4x 3 y 3 - 7x 3 y 2 - 7x 2 y 3 - Ux 3 y + Wx 2 y 2 - Uxy 3 + 13x 3 
+35x 2 y + 35xy 2 + 13y 3 - 30x 2 + Wxy - 30y 2 - 23x- 23y + 36 , 

Q(x, y) = {x- 1) (x + 1) (x 2 - 3x + l) (y - 1) (y + 1) (y 2 - 3y + l) . 



where 



and 



Tiling Crosses 

In how many ways can we tile a cross whose center is a 2 x 2 square and each of the four arms 
have length n? The answer is obtained by typing, in RITSUF: SeqCrossCsq(2,2,20,t) ;, getting 
(in 0.072 seconds!) that the number is 2B 2 (n) 2 , where 

£ gi( " r= ( t + l)(^-3 t+ l) ' 

n=0 



Incidentally, {S 2 («)} is http://oeis.org/A001654, the "Golden rectangle number" F n F n+ i, so 
the number of tilings of this cross is 2F 2 F 2 +1 , and it is a fairly simple exercise for humans to prove 
this fact. But we doubt that any human can derive, by hand, the answer to the analogous question 
for the cross whose center is a 4 x 4 square. The answer turns out to be B^n) 2 , where 

^ 3-t-5t 2 + 13t 3 -lH 4 -2t 5 + 2t 6 

^ ^ n)t ~ ~ (t - 1) (t 4 - 11 i 3 + 25 1 2 - 11 1 + 1) (i 4 + 7t 3 + 13 1 2 + It + 1) ' 

Even SBE needed 4.528 seconds to derive this formula after DZ typed: 



SeqCrossCsq(4,4,30,t) ; 



Once again, the approach is purely empirical. The computer finds all the tilings that completely 
cover the central square (or rectangle), possibly (and usually) with some tiles extending beyond. 
For each such scenario we have the state for each of the four arms of the cross. Then we use 
the previously computed (and remembered!) TM(a) matrices, find the corresponding entry in 
TM(a) n , and multiply all these four numbers. Finally we (meaning our computers) add up all 
these scenarios, getting the desired number. See the source code of SeqCross(a, b, N). 

Monomer-Dimer Tilings 

The beauty of programming is that, once we have finished writing a program, it is easy to modify it 
in order to solve more general, or analogous, problems. By typing ezraMDO ; in RITSUF the readers 
can find the list of analogous procedures for monomer-dimer tilings, where one can use either a 
1 x 2 a 2 x 1 or 1 x 1 tile, or equivalently, tiling with dominoes where one is not required to cover 
all the cells. 

For example, if A(n) is the number of ways of tiling the Knuth region (obtained by removing the 
central n 2 squares from an (n + 4) x (n + 4) array of squares) by dimers and monomers, the answer 
is much messier. It took SBE 15 seconds to discover that 



P(t) = -94i 30 + 1361 t 29 + 43975 1 28 -494267T 7 - 5787443 t 26 + 61 186056 t 25 + 266911158t 24 

-3200500450 t 23 -3505671568 t 22 +74767156291 1 21 -29007687275 1 20 - 796609853769 1 19 +823823428983 t 18 
+3924729557742 1 17 -4977782472712 1 16 -9040256915004 1 15 +11643454084810 i 14 +9751493606823 t 13 
-11693567793807 i 12 -4837640809485 1 11 +5123918478955 1 10 +1059903067708 t 9 -944330286322 t 8 -87120095554 1' 
+67451928324 i 6 + 657867045 1 5 - 1679236205 i 4 + 73176689 t 3 + 6962033 1 2 - 226706 1 - 10012 , 




where 



and 



Q(t) = (t - 1) (t 3 - It 2 + lit - 1) (t 3 + 7t 2 - 33t - 1) (t 3 - lit 2 + 7t - 1) • 

(t 3 + t 2 -3t-l) (t 3 -27t 2 + 107t -1) (t 3 + 3t 2 -t- 1) • 
(t 6 + 20t 5 + 55t 4 - 304t 3 - 337 1 2 + 8t + l) (t 6 - 37 1 4 - 76t 3 - 37 1 2 + l) 
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Weighted Counting 

What if instead of straight counting one wants to do weighted counting?. The weight of a domino 
tiling is defined to be fi# HorizontalTtles v # VertlcalTlles (i n the case of domino tilings) and 

h #HorizontalTiles v #VerticalTiles m #Mono m erTiles ^ Q thg cage Q f monomer-dimer tiling). 

For this we have the Maple package RITSUFwt, freely downloadable from 

http://www.math.rutgers.edu/~zeilberg/tokhniot/RITSUFwt, that also contains all the procedures of RIT- 

SUF. To get a list of the procedures for weighted counting, type: 
ezraWt () ; 

For example, to get the weighted analog of the Knuth problem, type: 
SeqFrameCwt(2,2,2,2,30,t,h,v) ; 
To see the output, go to: 

http : / / www . math . rutgers . edu/~zeilberg/tokhniot/ oRITSUFwt 1 , 

where one can also see statistical analysis. 
The Maple package ARGF 

The Maple package ARGF (short for "Analysis of Rational Generating Functions") downloadable 
from: 

http: //www. math. rutgers . edu/~zeilberg/tokhniot/ARGF , 

does automatic statistical analysis of random variables whose weight-enumerators are given by 
rational functions, like the one outputted by RITSUFwt, and whose procedures are also included in 
the latter. We urge the readers to look at 

http : / / www . math . rutgers . edu/~zeilberg/tokhniot/ oRITSUFwt4 , 

and the other output files in the webpage of this article 

http : / / www . math . rutgers . edu/~zeilberg/mamarim/mamar imhtml/r it suf . html , 

for examples. Here we use the methodology of [Z3][Z4]. 
Conclusion 

The deductive method in mathematics reigned for the last 2500 years. It is time to replace it by the 
inductive method. Often the inductive method is fully rigorous (like here), but other times, it may 
not be, but who cares?. Only trivial results can be proved fully rigorously, so let's not tie ourselves 
with these antiquated shackles, and explore mathematics experimentally! 
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